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Abstract o  Using  a  3-dimensional  formulation  for  an  aircraft’s 
dynamics,  the  required  controls  for  a  minimum  time- to- turn  are 

calculated.  Three  controls  are  usedi  (1)  angle  of  attack,  (2) 

♦ 

bank  angle,  and  (3)  thrust.  The  minimum  time-to-turn  solutions 
are  subject  to  varying  terminal  conditions  on  both  flight  path 
angle  and  heading  angle. 

In  general,  the  times  for  the  turns  are  not  greatly 
changed  by  varying  thrust/weight  ratios  or  the  final  flight 
path  angle.  Significant  effects  on  the  change  in  total  energy, 

final  altitude,  final  velocity  and  control  histories  are  noted 
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for  variations  of  the  above  parameters. 

Solutions  to  the  above  problem  are  accomplished  through 
the  use  of  Miele's  Sequential  Gradient-Restoration  Algorithm. 

1.  Introduction 

In  the  past  few  years,  optimization  techniques  have  been 
used  to  determine  turning  performance  of  high  speed  aircraft. 

Refs.  1-3  discuss  minimum-fuel  and  minimum-time  turns  at  con¬ 
stant  altitude.  Bryson  and  Hedrick  (Refs.  4-5)  have  extended 
the  analysis  of  these  types  of  turns  to  three  dimensions  by 
using  the  energy-state  approximation.  The  energy-state 
approximation  was  also  used  by  Beebee  (Ref.  6)  to  determine 
three-dimensional,  minimum-time  turns  for  p  hypersonic  rocket 
powered  aircraft.  Kelly  and  Edelbaum  (Ref.  7)  suggested  the 
use  of  asymptotic  expansions  in  conduction  with  energy-climbs 
and  energy- turns. 

Energy  approximations  are  very  helpful  in  studying  optimal 
aircraft  flight  since  the  order  of  the  state  differential  system 
is  reduced.  In  the  case  of  minimum- time-to-climb,  Bryson,  Desai,  & 
Hoffman  (Ref. 8)  have  shown  that  energy-approximations  give  good 
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results  when  compared  to  the  more  complex  dynamic  models.  Yet, 


there  is  a  need  to  check  the  accuracy  of  these  approximations 
when  used  to  determine  three-dimensional  tunning  performance. 

Preyss,  Willes,  Humphreys,  and  Roberts  (Ref.  9)  applied 
the  Pontryagin  Minimum  Principle  to  the  problem.  This  approach 
to  minimum  time  turns  results  in  the  usual  two-point  boundary 
value  problem  with  its  inherent  solution  difficulties.  The 
approach  used  in  this  paper  is  based  on  the  Sequential  Gradient- 
Restoration  Algorithm  developed  by  Miele  ,  Pritchard  and 
Damoulakis  in  Reference  10. 

2.  List  of  Symbols 

Xj  non-dimensional  distance  (x-direction) 

x2  non-dimensional  distance  (y-direction) 

x^  non-dimensional  distance  (z-direction) 

x^  non-dimensional  speed 

x^  flight-path  angle  (y) 

xg  heading  angle 

u^  pseudo  angle-of-attack  controller 

u2  bank  angle  controller 

u^  pseudo  thrust  controller 

L  non-dimensional  lift 

D  non-dimensional  drag 
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subject  to  the  differential  constraints 


jff=  *>(x.u,.)  (2) 

and  subject  to  the  boundary  conditions 
x(0)  -  x,  =  0 

(3) 

#<x#e)  =  0 

#f 

where  the  state  variable  x  is  an  n-vector,  the  control  variable 
u  is  an  ra-vector,  e  is  the  independent  time  variable,  is 
the  final  time  and  ^  is  a  q-vector  function  defining  the 
terminal  boundary. 

3.1  Equivalent  Fixed  Final  Time  Problem  The  Sequential 

Gradient-Restoration  Algorithm  was  formulated  for  the  final 
time  fixed.  This  transformation  is  always  possible  through  the 

use  of  a  parameter  (see  Ref.  10).  In  this  problem,  if  the 

time  e  is  redefined  by  the  transformation 

e  =  Tt  (4) 

where  T=e^,  then  equations  1-3  may  be  written  as^ 

I=t  (5) 


^Differentiation  by  the  independent  variable  t  is  denoted  by 
a  dot  ( . ) . 
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*  =  <p(x,U,Ttt) 


x(0)  -  x0  =  0 


»Mx,r)1  =  0 


(6) 


(?) 


(8) 


3.2  Aircraft  Equations  of  Motion  If  one  assumes  (a)  that  the 
aircraft  flies  a  coordinated  turn,  (b)  that  negligible  fuel  is 
consumed  during  the  maneuver,  and  (c)  that  the  error  due  to  the 
thrust  vector  rot  being  colinear  with  the  velocity  vector  is 
negligible,  then  the  equations  of  motion  in  the  velocity  axis 

'  I 

coordinate  system  are 

I 

^  =  rx^cos(x^)  COs(x6) 

*2  =  rx4C0S(x^)sin(x^) 


*2  =  -  r  x^sir.(x^) 


(9) 


*4  =  r(T-D-Wsin(x5)) 

I 

=  t(L cosCugJ-WcosCx^JJ/x^ 

=  rLsin(u2)/(x4cos(x^))  , 

The  initial  conditions  of  the  state  variable  x  are  assumed 

I 

given.  The  vector  function  ip  for  the  problem  to  be  studied 
is  given  by 
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I 


K'tx.r)  = 


,  |*5<1>  -  xf5 


!6(1)  - 


where  the  components  xj  through  x4  are  free  at  the  final  tine. 


The  equations  defining  L.D,’ and 


a  are 


L  = 

i 

D  =  (k1+k2a2)<Tx2 


!  a  =  e^3x3‘ 

! 

3.3  Control  Constraints  In  order  to  obtain  realistic 


results, 


>  it  is  necessary  to  keep  the  angle-of-attack.  the  thrust,  and 
the  "g”  loading  constrained  within  specified  limits.  These 
constraints,  are  specified  in  this  problem  by 

.  i 

0  *  a  £  1 

i 

°  4  “H  4  1  ;  (12, 

0  -  T  ^  T  !  ; 

1  max 

All  three  of  these  constraints  result  in  control  variable 

i  • 

inequality  constraints  and  may  be  handled  by  a  simple  sub- 

}  ; 

stitution  transformation.  For, equation  12(c)  the  transformed 


T;  =  winS> 


Equations  12(a)  and  12(b)  both  result  in  a  constraint  on  the 
angle-of-attack  •  These  two  constraints  can  be  seen  pictorially 
in  Figure  1# 

a 

x4 

Figure  1.  Angle-of- Attack  vs  Velocity  Constraint 

Fcr  velocities  below  the  comer  velocity, a”^,  the  transformation 
a  =  sin2(uj )  (14) 

is  used,  and  for  velocities  above  the  corner  velocity  the 
transformation 

L  =  sin2(u1)  (15) 

is  used.  The  same  control  u^  is  used  for  both  inequality 
conditions,  equations  12(a)  and  12(b),  since  these  two  conditions 
do  not  hold  simultaneously. 
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4.  Solution  Procedure 


The  solution  procedure  consists  of  alternate  applications 
of  a  gradient  phase  and  restoration  phases.  In  the  gradient 
phase,  nominal  functions  for  x(t),  u(t),  and  t  which  satisfy 
the  differential  constraints  and  boundary  conditions  are  assumed. 
By  minimizing  the  first-order  change  to  the  performance  index, 
variations  dx(t),  4u(t),  and  At  are  determined.  These  var¬ 
iations  are  added  to  the  nominal  functions  to  obtain  varied 
functions  x(t),  u(t),  and  r  which  produce  a  decrease  in  the 
performance  index.  Since  the  constraints  are  satisfied  only  to 
first  order  in  this  phase,  the  varied  functions  normally  do  not 
satisfy  the  non-linear  differential  equations  or  the  boundary 
conditions.  Consequently,  a  restoration  ohase  is  used  to  change 
x(t),  u(t),  and  f  so  that  the  non-linear  differential  equations 
and  boundary  conditions  are  satisfied. 

To  begin  the  restoration  phase,  the  varied  functions  x(t), 
u(t),  and  r  are  treated  as  nominal  functions.  By  minimizing 
the  least-square  change  of  u(t)  and  r  subject  to  the  quasi- 
linearized  differential  equations  and  quasil inear ized  boundary 

9 


m 


conditions,  a  new  set  of  variations  3x(t)t  i),  and  dr  is 
obtained.  These  variations  are  added  to  the  nominal  functions 
x(t),  u(t),  and  r  to  obtain  new  varied  functions  x(t),  ju(t), 
and  t  .  Tliis  step  is  performed  iteratively  until  the  differ¬ 
ential  equations  and  boundary  conditions  are  restored  to  some 
preselected  degree  of  accuracy.  The  cumulative  errors  in  the 
constraints  and  optimum  conditions  are  evaluated  using  the 
functionals^ 


-/! 

-/■ 

•'n 


N(i  -  <p  )dt  + 


N( 


(16) 


0  , 

+  N( 


NGC  +  O)  dt 

Jv 


*/ 


N(W)dt 


(f_-<0r)dt+(l*Hfcu)1 )  +  NU+g/u),  (17) 


Thus  P  must  be  reduced  to  an  acceptably  small  value  in  each 
cycle  of  the  Sequential  Gradient-Restoration  Algorithm.  When 
both  P  and  Q  are  acceptably  small  values,  convergence  has  been 


obtained. 


6  T 

The  norm  of  a  vector  z  is  defined  to  be  N(z)  =  z  z 


Initially,  nominal  solutions  which  satisfy  all  the  con¬ 
straint  conditions,  equations  6-8,  are  not  normally  known. 
However,  this  presents  no  problem,  since  only  the  restoration 


phase  need  be  applied  to  some  initial  guess.  In  this  way. 


feasible  nominal  solutions  can  be  generated  and  the  Sequential 
Gradient-Res toraticn  Algorithm  can  be  used  to  seek  the 


optimal  solution. 


5.  Necessary  Conditions 


As  in  other  approaches  to  optimizaion  problems,  this  method 


also  has  necessary  conditions  which  must  be  satisfied.  Both 


the  gradient  phase  and  the  restoration  phase  have  a  set  of 


necessary  conditions  which  result  from  application  of  opti¬ 


mization  criteria  such  as  the  Euler  equations  and  transversal ity 


conditions.  However,  the  particular  form  of  these  conditions 


is  a  consequence  of  the  numerical  technique  to  be  used,  A 


detailed  development  of  the  results  presented  in  the  next  few 


paragraphs  can  be  found  in  Refs.  10  -  11*  The  necessary 


conditions  are  as  follows i 


5.1  Gradient  Phas.? 


•  m  m  fp 

A  =  <p\k  +  <P^B  +  <P*  C 


*  -  -  ’’x* 


B*-«V 


*-A 


Adt  -  1 


where  A=4x/a,  B=du/at  C=dr/a,  and  a  is  the  step  size  for  the 


gradient  phase.  Equation  (18)  is  a  first-order  approximation 


of  the  differential  constraints  and  equations  (19)  -  (21)  are 


the  necessary  conditions.  The  boundary  conditions  are 


A(0)  =  0 


(  4  ^C)1 
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The  equivalence  to  the  transversal ity  condition  is 


(  A.  4  ^x«)1  =  0  (24) 

5.2  Restoration  Phase 

X  =  <p\  +  (p'h  +  <p^c  -  ( (p-  *)  (25) 

X  =  -4>xX  (26) 

B  =  -<Puk  (2?) 

C  =  -r^Xdt  (28) 

•'0 


where  A=dx/b,  B=4u/b,  C~Ar/\it  and  b  is  the  stepsize  for  the 
restoration  phase.  Equation  (25)  is  a  quasilinear  approximation 

of  the  differential  constraints  and  equations  (26)  -  (28)  are 
the  necessary  conditions.  The  boundary  conditions  are 


A(0)  =  0 

(29) 

(#  +  tf^A  +  =  0 

(30) 

The  equivalence  to  t^e  transversal ity  condition  is 

(A  +  4^xu)1  =0  (31) 
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The  matrices  <p£,  <p^f  and  <p£  are  given  by  the  same  expression  in 


both  the  gradient  and  restoration  phases.  For  the  matrix 
all  the  elements  ^(ij)  are  zero  except  for  the  following! 

(1.4)  =  rcos(x^)cos(xg) 

(1.5)  =  -rx^sin(x^)cos(xg) 

(1.6)  =  -rx^cos(x^)si n(xg) 

(2.4)  =  cos(x^)sin(x^) 

(2.5)  =  -rx^sin(x^)sin(x6) 

(2.6)  =  rx^cos(x^)cos(xg) 

(3.4)  =  -rsin(x^) 

(3.5)  =  -txj^cos(x^) 

(4.3)  =  -tD 

3 

(4.4)  =  -td 

x4 

(4.5)  =  -tWcos(x^)  (32) 

(5.3)  "  tLx^cos(u2)/x4 

(5.4)  =  -t( Lcos ( u2 ) -Wcos ( x$ ) ) /x J  +  tLx^cos(u2 )/x^ 

(5.5)  =  rWsin(x^)/x^ 

(6.3)  =  tLx  sin( u2 )/ x^cos ( x^ ) 

^  2 

(6.4)  =  -TLsin(u2)/x^cos(x^)  +  tLx  sin(u2)/x^cos(x^) 

2 

(6.5)  =  rLsin(u2)sin(x^)/x^cos  (x^) 

T 

For  the  matrix  <p  ,  all  the  elements  are  zero  except  for  the 
following! 


(4,1)  =  -tD. 


u, 


(4,3)  =  2rTmaxsin(u^)cos(u^) 
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(5*1)  =  tcos(u2)L^/x4 
(5,2)  =  -TLeinfugJ/x^ 

(6.1)  =  Tsin(u2)L  /x4cos(x^) 

A 

(6.2)  =  rLcos(t;  )/x4cos(x^) 

For  the  matrix^,  the  elements  are* 

(1.1)  =  x4cos(x^)cos(xg) 

(2.1)  =  x4cos(x5)sin(xg) 

(3.1)  =  -x4sin(x5) 

(4.1)  =  T-D-Wsin(x^) 

(5.D  =  ( Lcos  ( u2 )  -Wcos  ( x5 )  )/x4 

(6.1)  =  Lsin(u2)/ (x4cos(x^)) 

where,  for  x4  -a"^, 

L  =  k_,L  ;  L  =  2L/x4  ;  D  =  k3D  j  D  =  2D/x4 
Xj  J  x4  *3  ■>  *4 

L  *  k2  4sln(2ul’  '  Dut  =  2k2L1u/  4  °5) 

and  for  x4  -<r"^, 

2  ?  *  2. 

L  =  0  ;  L  =  0  ;  D  =  k-U,  x.-k2L  /  x4 ) 

Xj  .  x4  x3  j  j-  *+ 

D  =  2(kx  x4- k2L2/  x4)/x4  i  =  sin(2Ul)  (36) 

\  ‘  21t2llu/  4 


(33) 


(34) 
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6,  Numerical  Results 


T*e  numerical  results  presented  are  for  an  airplane  with 
the  following  nominal  characteristics! 

W  =  12,150  lbs.  kx  =  0.02 

"max  *  *2  ’  °-°5 

S  =  217  sq.ft.  C.  •  1.0 

rnax 

The  aircraft  was  assumed  to  be  flying  straight  and  horizontal 

immediately  before  entering  the  turn.  Thus,  the  flight  path 

\ 

angle  and  heading  angle  were  given  initial  values  of  zero. 

The  position  coordinates  Xj  and  x2  were  assigned  values  of  zero 

7 

and  the  altitude  was  specified  as  13.390  ft  . 

In  order  to  obtain  information  on  a  variety  of  turning 

situations,  the  thrust/weight  ratio,  initial  velocity,  final 

heading,  and  final  flight  path  angle  were  all  varied  as  shown 

in  Table  1.  Emphasis  was  placed  on  a  thrust/weight  ratio  of 

one  since  advanced  aircraft  are  expected  to  have  thrust/weight 

ratios  in  this  range.  The  initial  velocivy  was  varied  to  give 

values  approximately  130  fps  above  and  below  the  corner  velocity 

?Note  that  in  the  system  of  equations,  altitude  is  negative 
upward.  Therefore,  the  initial  altitude  in  the  computer 
simulation  was  equivalent  to  -13.390  ft. 
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Data  Set 

T/Vf 

X4(0) 

*6(t> 

Xlj 

;(T) 

1 

0.38 

621 

180 

( 

) 

2 

0.50 

3 

0.75 

4 

1.00 

5 

1.25 

6 

1.50 

1 

1 

1 

' 

7 

0.38 

903 

180 

0 

8 

0.50 

9 

o.75 

10 

1.00 

11 

1.25 

12 

1.50 

> 

t 

f 

t 

13 

1.00 

621 

180 

30 

14 

15 

15 

0 

16 

-15 

17 

f 

% 

f 

t 

-30 

18 

l.CO 

903 

180 

30 

19 

I 

15 

20 

0 

21 

-15 

22 

% 

f 

' 

\ 

f 

-30 

23 

1 .00 

621 

90 

30 

24 

15 

25 

0 

26 

-15 

27 

t 

> 

f 

> 

f 

-30 

28 

1.00 

903 

90 

30 

29 

15 

30 

0 

31 

-15 

32 

' 

f 

> 

’ 

' 

-30 

Table  1 .  Computer  Schedule 
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at  the  initial  altitude.  An  envelope  of  60°  for  the  final 
flight  path  angle  was  obtained  by  varying  this  final  condition 
from  -30°  to  +30°. 

In  addition  to  the  variables  mentioned  in  the  problem 
formulation,  two  others  which  are  important  in  combat  maneuvers 
were  evaluated  in  the  numerical  process.  The  change  in 
specific  energy  was  calculated  using 

dhe  =he(r)  -  hg(0)  (3?) 

where  he=-xy*-x£/2g.  Also,  a  turning  radius  for  the  130  turns 
was  calculated  using 

r=  ij^(x1(T)-x1(0))2+(x2(T)-x2(0))2+(x^(T)-x^(0))2J^  (38) 
6.1  Effect  of  Thrust/Weight  Ratio  To  determine  the  effect  of 
the  thrust/weight  ratio,  the  final  flight  path  angle  and  final 
heading  angle  were  given  values  of  0°  and  180°  respectively. 

Two  cases  were  consideredi  (1)  initial  velocity  below  the 
comer  velocity  and  (2)  initial  velocity  above  the  comer 
velocity.  In  both  cases,  the  thrust/weight  ratio  was  varied 
from  0,38  to  1.5»  These  two  cases  are  represented  by  data  sets 
1  through  12  in  Table  1.  The  corresponding  values  for  the 
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r'ft) 


Altf(ft)  Velf(fts)  ef (sec)  he(ft) 


1 

9,982 

757 

-496 

2356 

2 

10,452 

775 

10.8 

394 

2397 

3 

11,856 

796 

10.9 

2310 

2464 

4 

13,099 

828 

11.0 

4359 

2543 

5 

14,470 

859 

11.1 

6529 

2562 

6 

12,300 

794 

10.5 

2720 

2445 

7 

15,460 

784 

11.4 

-1029 

3053 

8 

16,212 

801 

11.5 

139 

3125 

9 

15,871 

809 

11.1 

-3 

2890 

1C 

16,158 

828 

11.1 

769 

2853 

11 

16,856 

855 

11.1 

2577 

3058 

12 

17,634 

886 

11.2 

3798 

2883 

13 

14,499 

768 

11.1 

4275 

2612 

14 

14,146 

782 

11.0 

4273 

2562 

15 

13.099 

828 

11.0 

4359 

2542 

16 

8,829 

791 

9.5 

-780 

2296 

17 

9,041 

805 

8.5 

-295 

2211 

18 

18,970 

773 

9.9 

2196 

3003 

19 

18,940 

788 

10.9 

2534 

3034 

20 

17,250 

855 

11.1 

2577 

2853 

21 

13.303 

811 

10.9 

-2515 

2779 

22 

6.510 

1035 

10.2 

-2910 

3495 

23 

14,491 

695 

6.4 

2603 

24 

13.792 

719 

6.0 

2434 

25 

13.U9 

741 

5.7  ' 

2275 

26 

12,571 

759 

5.5 

2144 

27 

12,107 

776 

5.3 

2063 

28 

15.451 

945 

6.5 

3264 

29 

14,588 

973 

6.5 

3274 

30 

13,600 

1006 

6.4 

3280 

31 

12,519 

1040 

6.3 

3277 

32 

11,462 

1072 

6.2 

3264 

Table  2,  Computer  Results 


1? 


! 


t 

change  in  specific  energy,  turn  radius,  and  the  final  values 
for  altitude,  velocity  and  time  are  given  in  Tablp  2.  As*  can 

i 

be  seen  from  Table  2,  the  thrust/weight  ratio  has  little 
effect  upon  the  time  to  turn.  However,  it  does  affect  final 

:  i 

velocity,  final  altitude,  change  in  specific  energy,  'and  turn 
radius. 

When  entering  the  turn  with  velocity  less  than  the  corner 
velocity  (data  sets  1  -  6),  an  increase  in  thrust/weight  ratio 

1  ’  i 

produces  an  increase  in  final  altitude,  final  velocity,  and 

specific  energy.  However,  the  results  for  a  thrust/weight  ratio 

j 

of  1.5  are  not  consistent  with  these  general  trends.  This 

variation  from  the  trend  is  a  result  of  throttling.  Until  the 

thrust/weight  ratio  reaches  1.5»  turns  are  made  a  T=T  .  For 

,  max 

a  thrust/weight  ra+io  of  1.5*  sufficient  acceleration  exists 

» 

for  the  aircraft  to  fly  to  the  comer  velocity,  reduce  the 

thrust  to  zero  to  remain  at  the  corner  velocity,  and  finally 

.  *  . 

re-ayply  full  thrust  to  accelerate  on  the  n  boundary  to 

max  i 

complete  the  maneuver.  The  thrust  programs  for  this  caae  are 
shown  in  Figure  2.  1 
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When  entering  the  turn  at  a  velocity  greater  than  the 

comer  velocity  (data  sets  7  -  12),  an  increase  in  thrust/ 

•  »  1 
i 

weight  ratio  produce?  an  increase  in  final  altitude,  final 

*  -  * 

velocity  and,  in  most  instances,  specific  energy.  The  vari- 

,  I  ; 

ations  from  this  trend  are  also  caused  by  throttling.'  As  shown 
in  Figure  3*  throttling  occurs  when  the  the  thrust/weight  ratio 

;  i  ‘ 

is  0.75  or  greater.. 

' 

i 

In  the1 two' cases  considered,  the  aircraft  must  be  flown 

i  i 

! 

J 

along  the  stall  and/or  the  n  boundery.  Velocity  change*, 
i  •  max  (■ 

!  , 

I 

are  made  in  the  early  part  of  the  turn  to  drive  the  aircraft 

■  '  '  ’  '  ' 

velocity  toward  the  comer  velocity.  If  the  turn  is  entered 

t  • 

at  a  velocity  less  than  vc,  the  aircraft  is  flown  along  the 

1  i  . 

i 

stall  boundary  until  the  corner  is  reached.  The  remainder  of 

I 

the  turn  is  made  on  the  ^max  boundary.  When  entering  a  turn 

\ 

at  velocity  greater  than  vc,  the  aircraft  is  flown  along  the 

n:  boundary  until  the  comer  velocity  is  reached.  The 

max  i  i  i 

aircraft  remains  on  or  about  the  corner  velocity  for  a  short 

■  •  i 

time  and  then  completes  the  maneuver  by  accelerating 
on  the  n  boiindary.  In  every  case,  Maximum  throttle  is 
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applied  prior  to  reaching  the  comer  velocity.  As  the  thrust/ 
weight  ratio  increases,  the  maximum  throttle  is  initiated  when 
closer  to  the  comer  velocity.  Whenever  throttling  occurs,  the 
aircraft  can  remain  at  the  comer  velocity  for  a  longer  period, 
and  consequently,  the  time  to  turn  is  reduced. 

The  necessary  bank  angles  for  these  two  cases  are  shown  in 
Figures  4  and  5«  As  can  be  seen  from  these  figures,  minimum 
time  turns  have  a  given  bank  angle  to  start  the  turn  and  a 
continuously  varying  bank  angle  throughout  the  turn.  Also, the 
bank  angle  is  not  zero  at  the  end  of  the  turn.  Since 
bank  angle  is  a  control  variable,  it  can,  by  assumption, 
attain  ai.y  value  instantly j 

6,2  Effect  of  Final  Flight  Path  Angle  To  determine  the  effect 

of  the  final  flight  path  aj.^le,  the  thrust/weight  ratio  was 

given  a  value  of  1.0.  Both  180°  and  90°  turns  with  entry  speeds 

above  and  below  v  were  considered.  These  cases  are  represented 

c 

by  data  sets  13  -  32  in  Tables  1  and  2, 

6,2,1  180°  Turns  As  shown  in  Table  2  (data  sets  13  -  22),  the 
final  flight  path  angle  generally  has  little  effect  upon  the 
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time  to  turn  180°.  However,  final  altitude,  final  velocity, 
turn  radius,  and  change  in  specific  energy  are  strongly  affected. 
When  entering  the  turn  with  a  velocity  less  than  v  ,  throttling 

v 

occurs  for  the  negative  final  flight  path  angles  (Figure  6), 

In  all  cases,  however,  the  turn  is  started  at  T=T  •  When 

max 

entering  the  turn  with  velocity  greater  than  v  ,  throttling 

V 

also  occurs.  In  this  case,  the  turn  is  started  with  T=0  and 
power  is  added  after  a  period  of  time  (Figure  7). 

Specifying  the  final  flight  path  angle  has  a  definite 
effect  upon  the  type  of  maneuver  required.  For  example,  if  the 
aircraft  is  flying  at  a  velocity  less  than  v  ,  and  the  final 

V 

flight  path  angle  is  required  to  be  positive,  the  aircraft  will 
start  into  a  turning  dive  to  gain  velocity  and  end  in  a  climbing 
turn.  An  opposite  maneuver  is  required  when  entering  the  turn 
at  a  velocity  greater  than  v  and  when  a  negative  final  flight 

V 

path  angle  is  specified.  This  would  require  a  climbing  turn 
through  the  first  part  of  the  turn  and  a  turning  dive  at  the 
end.  The  bank  angles  required  for  these  maneuvers  are  shown  in 
Figures  8  and  9« 
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6,2,2  90°  Turns  For  90°  turns,  the  time  to  turn  decreases  as 
the  final  flight  path  angle  decreases.  This  trend  holds  true 
for  initial  velocities  above  and  below  the  comer  velocity 
(data  sets  23-32).  As  the  final  flight  path  angle  decreases, 
the  aircraft  takes  advantage  of  the  acceleration  of  gravity. 

Thus,  in  the  shorter  duration  turns,  load  factor  is  more 
important  than  aircraft  velocity.  There  is  no  exchange  of 
altitude  for  velocity  to  achieve  an  aircraft  velocity  near  v  • 

V 

Instead,  the  aircraft  is  flown  on  the  n  boundary  with  the 

max 

bank  angles  shown  in  Figures  10  and  11. 

6.3  Comment  on  the  Sensitivity  of  Control  Variables  Analysis 
of  the  results  of  the  digital  simulation  shows  that  changes  in 
the  control  variables  may  have  relatively  little  effect  on  the 
minimum  time  to  turn.  Although  these  characteristics  are 
desirable  for  a  pilot  performing  such  maneuvers,  they  provide 
minor  numerical  difficulties  in  obtaining  a  solution. 

The  computer  simulations  show  that  the  minimum  time  to 
turn  is  not  very  sensitive  to  the  bank  angle  control.  Relatively 
large  changes  in  the  bank  angle  program  (on  the  order  of  several 


degrees  between  itterations)  caused  virtually  no  change  in  the 


time  turn.  Precise  control  of  the  aircraft  in  bank  angle, 
therefore,  is  not  required  to  achieve  good  results. 

The  throttling  of  the  thrust  is  shown  as  a  time  varying 
function  over  rather  short  time  periods  (on  the  order  of  a  few 
seconds).  As  the  thrust/weight  ratio  gets  large  (greater  than 
or  equal  to  one),  the  thrust  program  approaches  a  bang-bang 
solution.  For  intermediate  values  of  thrust/weight  ratio  (like 
0.75) »  the  thrust  must  be  applied  well  before  the  corner  velo¬ 
city  is  attained  due  to  the  lag  from  relatively  low  acceleration. 
Higher  thrust/weight  ratios,  therefore,  provide  for  more  easily 
applied  thrust  programs  which  is  true  of  bang-bang  solutions. 

As  shown  earlier,  the  angle-of-attack  is  always  limited  by 
the  maximum  load  limit  or  stall  angle-of-attack.  These  two 
boundaries  are  natural  ones  to  a  pilot  and  it  is  not  difficult 
for  a  pilot  to  fly  along  these  boundaries. 

6,4  Comments  on  the  Algorithm  A  characteristic  of  the  numerical 
algorithm  is  the  ease  with  which  to  choose  the  initial 
nominal  trajectory.  For  the  examples  used 
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in  this  paper,  three  constant  values  of  the  control  are  assumed 


and  used  to  numerically  integrate  the  non-linear  equations  of 
motion.  After  the  initial  restoration  was  completed,  each  inter¬ 
mediate  sub-optimal  solution  always  met  all  the  constraints  on 
the  problem.  Therefore,  all  solutions  are  feasible  and  can  be 
compared  directly. 

Because  of  the  insensitivity  of  the  optimal  solution  to  the 

—4 

controls,  a  very  strict  convergence  criterion  of  Q-10  was  not 

-2 

used.  Solutions  for  which  Q-10  displayed  very  small  changes 
in  the  cost  functional  I=r. 

As  a  final  check  on  the  correctness  of  the  optimal  solution, 
the  control  histories  derived  by  the  Sequential  Gradient-Restor¬ 
ation  Algorithm  were  used  to  numerically  integrate  the  non-linear 
differential  equations.  In  every  case,  the  solutions  obtained 
from  this  final  integration  process  matched  the  optimum  solution 
derived  from  the  Sequential  Gradient-Restoration  Algorithm, 

?.  Summary 

A  method  which  utilizes  the  generalized  equations  of  motion 
to  determine  three-dimensional  minimum-time  turns  for  aircraft 
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has  been  presented*  The  problem  was  formulated  for  specified 

final  heading  and  flight  path^igles.  The  cnalysis  was  restricted 

to  subsonic  flight  only  by  the  expression  used  for  drag.  The 

solution  procedure  and  computational  procedure  were  based  upon 

the  methods  developed  in  Reference  10  • 

The  system  of  equations  for  the  problem  were  programmed  in 

ALGOL  for  the  Burroughs  5500  computer.  The  computer  simulation 

was  used  to  obtain  information  on  a  variety  of  turning  situations* 

Initial  velocity,  thrust/weight  ratio,  final  heading,  and  final 

flight  path  angle  were  found  to  have  a  strong  influence  on  the 

required  maneuver.  Minimum  time  turns  are  always  made  on  the 

stall  and/or  n  boundaries.  During  the  turn,  throttling  from 
max 

T=T  to  T=0  or  the  reverse  is  frequently  required  for  thrust/ 
max 

weight  ratios  -  0.75*  For  low  thrust/weight  ratios  an  altitude- 
velocity  exchange  is  used  to  approach  the  corner  velocity. 
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Figure  2.  Percent  Thrust  for  Various  T/Vf 
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Figure  3«  Percent  Thrust  for  Various  T/» 
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Figure  5»  Bank  Angle  for  Various  T/Vf 
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Pigure  6.  Percent  Thrust  for  Various  Final  Flight  Path  Angles 
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Figure  ?.  Percent  Thrust  for  Various  Final  Flight  Path  Angles 
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Figure  9.  Bank  Angle  for  Various  Final  Flight  Path  Angles 
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